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We re-analyze the direct J/ip production processes at the Fermilab Tevatron in view of the recent observation at the 
B-factories, where both J/ip inclusive and exclusive production rates are found to be about an order of magnitude 
larger than the leading order estimates of non-relativistic QCD. The charm quark fragmentation to J/4', which is the 
dominant color-singlet process at high pT, is normalized by the B-factory measurements. The process receives further 
enhancement due to the charm sea contribution which has so far been ignored in most analyses. After summing up all 
sub-process contributions, we find that the color-singlet mechanism alone can account for 20% to 90% of the observed 
direct J ftp high pT production. The polarization rate of the directly produced J/i/i is sensitive to the fraction of the 
color-octet contribution, which is employed to fill the gap between the color-singlet prediction and experimental data. 
With a bigger ii'-factor for the charm quark fragmentation probability, we envisage a smaller matrix element for the 
color-octet "^5^**' state, and this can be examined at the LHC in near future. 



1. INTRODUCTION 

In the history of quarkoniuni physics, the so-called color-singlet model (CSM) was widely employed in the study of 
heavy quarkonium production and decays [jI . Although CSM made a great success in the description of quarkoniuni 
production and decay in the past, the Tevatron data 0,01 on large-pT J/ifi production collected in the 1992-1993 run 
severely challenged it, as the data differ very much from the leading order (LO) CSM predictions in both normalization 
and pt distribution. 

One of the developments in past decade is the so-called nonrelativistic quantum chromodynamics (NRQCD) 
which provides a plausible explanation for the CDF measurements on direct J/ -0 and i/j' production at the Fermilab 
Tevatron Within the framework of NRQCD, the surplus charmonia production in excess of the CSM prediction 
attributes to the color-octet contribution. Since the magnitudes of the universal color-octet matrix elements are 
determined by fitting to the experimental data, to find an independent color-octet signature is necessary to further 
confirm the applicability of NRQCD to charmonium production. The color-octet scenario encounters some difficulties 
while it is confronted with the new experimental data [H]- The most striking puzzle probably is the absence of 
transversely polarized J/ ijj and i/'' in experiment 0, 0] ■ 

Recently, BaBar and Belle data indicate that the experimental measurement 0, @, [l3l of the J/ ip inclusive and 
exclusive production rates at the B-factories greatly overshoot the leading-order theoretical calculations in CSM [lli| . 
Because the J/ips are found to be produced in accompany with charmed mesons, very large higher-order corrections 
to the charm quark fragmentation process should be inevitable, though fragmentation mechanism may not work well 
at or beyond the next-to- leading order (NLO). According to the analyses in Refs. 12, l3|, + e~ J/tp + c + c 



process requires a /iT-factor of 2.6l^'^ to 7.2 ± 2.9 to mimic the higher order corrections. It was found that with a 
constant ii'-factor of around 4 one can reproduce both the inclusive J/^ production momentum distribution and the 
large branching ratio of J/ip + c + c production observed at the B factories 0, [13, Indeed, very recently large 



NLO corrections to the exclusive double charmonia production [15| and the J/ip + c + c production processes [16| 
are obtained. Note that the employed ii'-factor in this work represents all corrections, rather than merely the NLO 
QCD one. 

In the following, by taking into account the charm-sea induced processes and a iiT-factor obtained in the B-factory 
in representing various higher-order corrections to the charm quark fragmentation mechanism, we re- visit J/^ surplus 
production issue at the Fermilab Tevatron. 
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Figure 1: Direct J/ij; production diflPerential cross-sections for various color-singlet processes versus the Run-I experimental 
data from the Fermilab Tevatron. 



2. CALCULATION AND RESULTS OF THE J/^ HIGH pt PRODUCTION 

Schematically, the quarkonium hadroproduction process, A -\- B ^ + X, can be formulated in a factorized 
form, 

dcj{A + B ^ H{pt) + X)=y2 [ dXa [ dXb fa/A{Xa,Q) 

„ I, Jo Jo 



X fb/B{xb,Q) da{a + b ^ H{pt) + X) , 



(1) 



where a and b are the incident partons within the colliding hadrons A and B, respectively; fa/ a and fb/B are parton 
distribution functions at the scale of Q^. In larg e transverse momentum, the charmonium hadroproduction are 
dominated by the fragmentation sub-processes [l7| , which can be expressed as 

da{a + H{pt) + ^) = / dz da{a + b^ c{pt/z) + X) Dc^h{z, Q) , 

. Jo 



(2) 



where c denotes the fragmenting parton, Dc^h{z,Q) is the parton fragmentation function. The evolution of the 
fragmentation function Dc^h{z, Q) in regard of scale Q in Eq. ^ is realized by solving the Altarelli-Parisi (AP). 

In our numerical calculation, J/?/; radial wave function at the origin i?(0) is obtained from its Icptonic decay width 
T{J/t/j e+e~) = 5.55 x 10^^ GeV [18], and the CTEQ6L parton distribution functions are employed To avoid 
the large logarithms from high order corrections, the fragmentation, factorization and rcnormalization scales are set 
to be the J/ip transverse mass, i.e. tot - 



Pt 



Figure [T] shows the transverse momentum distribution of the directly produced J/i/j at the Fermilab Tevatron, 
where eight main color-singlet sub-processes and their sum are presented separately. In the literature, the analytic 
results for most of the CS processes exist, and we find numerically agreements between our results and those in 
previous studies. Different from Refs.[2oi 21 1, here the gluon fragmentation process overshoots the charm quark 



fragmentation process, due to the reason that the induced charm quark fragmentation in the gluon fragmentation 
evolution equation is taken into account in our calculation. From the figure, one can easily find that two charm-sea 
induced sub-processes contribute dominantly to the direct J/?/' production, which in fact partly represent the next- 
to-next-to-leading(NNLO) order g-\-g^J/il)-\-X process. Although the importance of the second process at very 
high transverse momentum was mentioned in Ref. (2^ . the inclusion of the first one is necessary at moderate pt 
region. 

Enforcing a factor of 4 to the charm quark fragmentation function and an overall factor of 2 to the hard scattering 
processes, we obtain the re-normalized differential cross sections, as shown in Fig. [2ja). Here, the difference between 
experimental data and theoretical calculation in GSM reduces to about a factor of 4. If we assume that the remaining 
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Figure 2: NRQCD prediction for the direct J/^j) production at the Tevatron Run I energy, where the color-octet matrix element 
is adjusted in order to fill the gap between CS prediction and experimental data, (a) The normalization factor for the charm 
fragmentation process is 4 x 2, 4 for is'-factor and 2 for hard-scattering cross-section; (b) the normalization factor is taken to 
be 8 X 2. (c) The J/tj} polarization variable a versus transverse momentum. 



gap attributes to the color-octet contribution, the matrix element of < Ops'!* |0 > may diminish 1/4 from the initial 



fitted value 23|. In the case oi K — A: one cannot find a consistent explanation for both CDF and Belle results. 
Therefore, the larger X-factor of 7.2 ± 2.9, which is the upper limit accommodated by the Belle data is envisaged. 
In this case, the color-singlet mechanism alone explain the data in large transverse momentum region, as shown in 
the Fig. Hb). 

Figure [2jc) exhibits our calculation result on J /tp polarization distribution versus transverse momentum. Here, 
the higher excited state feed down effects are neglected. Of this figure, the vertical axis denotes the conventional 
measure of polarization, defined as: a = [ar ~ 2aL)/{crT + '^o'l)', the horizontal axis denotes the J/4> transverse 
momentum. The upper and the lower limits, a = 1 and a — ^l, correspond to the solely transverse and longitudinal 
polarization cases respectively, while a = corresponds to the unpolarized situation. In our result, because the quark 
fragmentation process gives unpolarized J/ip in the relativistic limit 2J] , the total polarization of directly produced 
J/ip in high px region tends to be unpolarized. This is completely different from the predictions in color-octet model 
and in previous calculations in CSM [2o| . 



3. CONCLUSIONS 

In summary, we have recalculated the direct J/tp high-p^^ production rate in color singlet mechanism and find 
that the charm sea induced processes contribute dominantly in the direct J/tp production while Pt> 4 GeV at the 



Fermilab Tevatron. In recently the next-to-leading order QCD correction to process (5) is realized l25|,|26|. However, 
as it is well-known that in J/ip high px hadroproduction the dominant processes are those beyond the NLO ones. 
So, only the full NLO calculation is not enough for the issue discussed in this work. 

Based on the iC-factor scenario we find that the Tevatron data can be almost accommodated by the color-singlet 
production mechanism. We assign if-factors of 4 and 8 respectively to the charm quark fragmentation process, and 
a factor of 2 to the parton-parton hard interaction processes. Quantitatively these enhancement factors may vary 
from process to process, but for a qualitative analysis we think a constant number is still meaningful. In the case of 
if-factor equals to 8, the CDF data can be simply understood in the color-singlet mechanism. 

We have also evaluated the polarization rate of directly produced charmonium in color-singlet mechanism. We 
find that at large transverse momentum, the J/tps will be statistically unpolarized. If this is confirmed by future 
experiment, then we may conclude that the gap between our calculation result and experimental data should be 
eliminated by giving a bigger X-factor to CS processes, but rather filled by the color-octet component. This is in 
fact somewhat favored by the Tevatron II preliminary measurement (27| . 
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